Spin Susceptibility of an Ultra-Low Density Two Dimensional Electron System 
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We determine the spin susceptibility in a two dimensional electron system in GaAs/AlGaAs over 
a wide range of low densities from 2xl0 9 cm -2 to 4xl0 10 cm -2 . Our data can be fitted to an 
equation that describes the density dependence as well as the polarization dependence of the spin 
susceptibility. It can account for the anomalous g-factors reported recently in GaAs electron and 
hole systems. The paramagnetic spin susceptibility increases with decreasing density as expected 
from theoretical calculations. 
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PACS numbers: 73.40.-c, 71.27+a, 71.30+h 

The low-density ground state of a degenerate elec- 
tron system is one of the oldest questions of many par- 
ticle physics. In two dimensions, it is expected that 
in the dilute limit, the electron Fermi liquid under- 
goes a phase transition to a solid known as the Wigner 
crystalQ, 0, Furthermore, driven by the exchange 
coupling between electrons, a ferromagnetic state may 
arise at densities slightly above the critical density of 
Wigner crystallization^ H, El 0- None of these phases 
have yet been detected due to the lack of sufficiently high- 
quality, low-density specimens. Recent measurements of 
relevant system parameters such as the spin suscepti- 
bility x an d the effective g-factor g* have shown con- 
siderable deviations from their standard, non-interacting 
valuesjl H [H [H E3 This is believed to be the re- 
sult of strong interactional 0, but disorder-led enhance- 
ment cannot be ruled out. Such explorations of the pre- 
transition regime provide valuable data against which to 
check the same theoretical calculations that determine 
the transition point to either a ferromagnetic state or 
to the Wigner solid. A possible connection between the 
ferromagnetic properties and the metal-to-insulator tran- 
sition (MIT) in low density 2DESs has added to the com- 
plexity of the problem and prompted intense recent stud- 
ies on this subject in Silicon MOSFETsHEiiElElIli. 
These results generally lend support to the possibility of a 
ferromagnetic state. On the other hand two groups have 
recently reported anomalous density dependencies of the 
g-factors in the GaAs/AlGaAs electron and hole systems 
that are at odds with results in Silicon MOSFETs and 
disfavor the predicted ferromagnetic transition [ill fl^| . 

In this paper, we report measurements of \ in a vari- 
able density 2DES of exceedingly high quality to unprece- 
dented low densities. In addition to the density depen- 
dence of Xj our measurements determine, for the first 
time, the explicit polarization dependence of x- This po- 
larization dependence can account for the " anomalous" 
g-factors recently re por ted in GaAs 2DESs and 2D hole 
systems (2DHG) [H H3 . 



Our specimen is a heterojunction-insulated gate field 
effect transistor (HIGFET). The specimen consists of a 
(001) GaAs substrate, overgrown by molecular beam epi- 
taxy with 0.5/xm of GaAs, followed by a 200-fold su- 
perlattice of lOnm of GaAs and 3nm of Alo.32Gao.68As. 
Subsequently, 2/im of GaAs are deposited as a channel, 
followed by 600nm of Alo.32Gao.6sAs as an effective in- 
sulator, and capped by a heavily doped GaAs n + layer, 
serving as a top gate. The specimen is processed into 
a 600 fim square mesa. Sixteen Ni-Ge-Au contact pads 
are spaced evenly along the edges of the mesa using stan- 
dard photolithography. One corner pad provides the con- 
tact to the top gate, which allows for a continuous and 
in situ change of the 2DES density. The density range 
available for measurement extends from 1.7xl0 9 cm 2 to 
6.4xl0 10 cm 2 . The interaction parameter r s , defined as 
the ratio of the inter-electron spacing to the Bohr radius, 
a/as=l/Y / 7mas, spans 2.2 < r s < 13.4. Screening of 
the e-e interaction via the top metallic gate can be ne- 
glected, since its distance exceeds the electron spacing 
by a factor of 4.4 to 27 in our experiments. The MIT 
occurs at 2xl0 9 cm 2 , which is the lowest transition den- 
sity ever reported in a 2D system, attesting to the low 
disorder of the specimen. Our FET gives us the unique 
advantage of wide density tunability within one specimen 
and its reproducible gate- voltage/density relation lets us 
sweep density at fixed magnetic field, which is an essen- 
tial method in our x determination. 

Our measurements were performed in a dilution refrig- 
erator equipped with a rotating sample platform reaching 
a base temperature of 30mK. A standard low frequency 
(3-23Hz) lock-in technique was used with excitation cur- 
rents ranging from lOOpA to lOOnA. All experiments were 
performed during one cooldown, facilitating quantitative 
comparison between our data. 

In a normal Fermi liquid, the spin susceptibility x 
=dAn/dB= g*/is / o/2, where g* is the effective g-factor, 
p is the density of states (DOS) at the Fermi level, and 
An=n^-n-k In a 2D system, p — m*/nh 2 , therefore 
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X = g*m* ^b/27t?i 2 . Generally, theories find that x in- 
creases with growing interaction due to spin-exchange 
coupling. To investigate the experimental situation, we 
measure the spin susceptibility x as a function of the 
2DES density n. We express it as a relative spin sus- 
ceptibility x/xo = m *9* / m b9b, where m(,=0.067m e and 
gh=0.44 are the band values of mass and g-factor in GaAs 
and xo the Pauli susceptibility determined by these band 
values. In the remainder of the paper, we use values for x 
and m*g* normalized to xo and nibgt, respectively, so that 
X= m*g*. We employ two different methods to measure 
m*g*. First, we follow and extend the tilted field method 
introduced by Fang and Stiles in Silicon MOSFET0. 
Secondly, we follow the parallel-field method recently uti- 
lized by Refs. 0,0] to derive m*g* from the full polar- 
ization condition of the 2DES. 

In a magnetic field, spin-up and spin-down electrons 
form two separate sets of landau levels. As the mag- 
netic field is tilted, the two sets of landau levels move 
with respect to each other. The energy diagrams are 
schematically shown as insets in Fig. Qfa). Solid and 
dotted lines represent spin-up and spin-down Landau 
levels, respectively. The spacing between Landau lev- 
els, hoj c = eKBperp/m*, depends on the perpendicular 
component of the field. The shift between both sets, on 
the other hand, is the Zeeman energy, AE Z = g*^B^tot, 
which depends on the total field. By adjusting the tilt 
angle, 9 and the total field, B tot , hcu c and AE Z can be in- 
dependently changed. Particularly useful configurations 
arise when g*^sBtot=i?iw c , where i is an integer or a half- 
integer. At half-integer configurations, the Landau levels 
from both spins interleave and form a set of uniformly 
spaced levels. At integer configurations the Landau lev- 
els coincide and form again a set of uniformly spaced 
levels, however, this time with double spacing compared 
to the half-integer case (see insets Fig. EI a)). 

The magneto-resistance of the 2DES reflects the con- 
figuration of the energy levels. Distinctive signatures 
from different configurations are observed when the 
2DES density is swept at fixed Utot and 9. B perp , and 
hence 9, is accurately determined from the period of the 
oscillations and the Landau level degeneracy eB perp /h. 
At half-integer configurations, the depths of successive 
minima are just equal. This is also the case for the 
integer configurations, but here each second minimum 
disappears. The different configurations are realized 
only transiently at particular densities no that satisfy 
g* [iBB to t=ifiu c , equivalent to m*g*(no)=ie?icos#/yLtB, 
since m*g* depends on density. Integer configurations 
are easily identified by the disappearance of a minimum. 
At half-integer configurations, the depth of neighboring 
minima interchange strength at no, which is identified by 
the crossing point of two smooth envelopes drawn along 
alternate minima, see Fig. da) for 85.71°. The trace 
at 88.40° of Fig. da) shows, as examples, several con- 
figurations as the density is swept and indicated by the 
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FIG. 1: (a) Density sweeps at B perp =0.07T and 9 = 85.71° 
(thin) and 88.40° (thick) respectively. Events correspond- 
ing to g*/^sB t0 i=i?iu; c; where i=l/2, 3/2, 2 and 5/2 are indi- 
cated in the plot with level diagrams nearby as insets. Dashed 
lines illustrate the procedure to determine no at half-integer 
configurations, (b) Linear dependence of coefficients of the 
power law fits of Fig. [5] on their index, i. (c) Magneto- 
resistances in parallel-field method for different densities in 
unit of 10 10 cm~ 2 . Positions of full polarization fields, B p , are 
indicated following Ref.[lj|. (d) B p as a function of density. 



neighboring diagrams. Increasing 9 slowly and tracking 
the indices carefully, we are able to identify and label 
events that belong to configurations with index 1/2, 1, 
3/2, 2, and 5/2. The product m*g* is calculated accord- 
ing to m*g*=ie?i cos 6>//is. The solid symbols in Fig. [3 
represent data derived with this method, using different 
symbols for different indices. 

Clearly, m*g* increases with decreasing density and 
increasing index i. More strikingly, for each fixed index 
i, m*g* displays a power law dependence on n within the 
measured range, indicated by the dashed lines on a log- 
log scale. This dependence is best developed for i= 1/2 
and 3/2. Furthermore, the parallel lines suggest a single 
exponent. We therefore fit a power law dependence to 
the data of i=l/2 and use the exponent as a constraint 
in fitting data of the other indices. The coefficients from 
these five fits display a very good linear dependence on 
i, as seen in Fig.^b). This gives us an empirical equa- 
tion, m*g* =(2.73+0.66*i)n -0 - 4 , where n is in unit of 
10 10 cm~ 2 . This equation describes all our data points 
remarkably well. 

Before discussing the implications of such an empir- 
ical equation, we proceed with the second method of 
determining m*g*. With increasing in-plane magnetic 
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field, B||, the polarization, P=(n^-n^)/n of the 2DES in- 
creases and saturates at unity at a threshold field B p 
when g*/isB p /2=E^. The derivative of An (=Pn) with 
respect to By is x( = m*g*). Assuming m*g* to be 
independent of P it follows that m*g*=27T?i 2 n//iBB p . 
As asserted by several groups, full polarization of the 
2DES is signaled by the onset of an exponential behav- 
ior in a parallel-field magneto-resistance experiment 0, 
UtI Hil fl3. |20| We have performed such experiments on 
our HIGFET for different densities (Fig. ^c)), deter- 
mined Bp according to Ref. 1 an d show their values in 
Fig. [2d). Within experimental error, B p is independent 
of the current direction relative to B||. Translated into 
m*g* the results are plotted as open circles in Fig. 
Unlike the solid data points from the first determination 
of m*g*, these new data exhibit a non-monotonic de- 
pendence on density and are consistently larger in value 
than the solid data points, with the discrepancy increas- 
ing with increasing density. 

The key to reconciling the discrepancy between these 
data sets is to recognize the dependence of m*g* 
on the polarization, P. In the tilted field method, 
a higher index i implies a higher Zeeman energy 
and therefore a higher degree of 2DES polarization. 
In the Fermi liquid limit (B perp =0), P=(n^-n^)/n = 
AE z p/2n={g*ii B B tot )(m*/2iTh 2 )/n. The tilted field 
method assumes that the introduction of a small B perp 
does not alter the values of g* and m*. Using the rela- 
tion g*/xsBt ot =i?kj c , a straightforward calculation yields 
P=ieB peT . p /nh. This direct relationship between P and 
index i establishes an explicit P-dependence of m*g*. 
Substituting P for i and using B perp =0.07T[2l], we ob- 
tain x=m*g*=(2.73+3.9Pn)n-°- 4 from our data. Most 
remarkably, for fixed density n, the susceptibility x m_ 
creases linearly with increasing polarization. This mono- 
tonic increase reflects an increasing spin-exchange energy 
with increasing population of like-spins. 

Quite clearly, the assumption of a P-independent x, 
assumed in Refs. to derive m*g* from B p , does 

not apply. Our empirical interpolation equation can 
be used to make contact between the data obtained 
from the tilted-field method and from the parallel field 
method. Their relationship is best discussed referring to 
the inset of Fig. [21 This diagram shows, as an exam- 
ple, the evolution of the net spin An=Pn=n^-n^ with 
B to t for a fixed density, n=lxl0 10 cm~ 2 . By definition, 
the slope of An(B) is the susceptibility x= m*g*. The 
parallel-field method determines An=n (=lxl0 10 cm~ 2 ) 
at B p =4.9T. The assumption of a B-independent m*g* 
implies a linear rise of An with B, shown as a straight line 
in the inset. One value, B p =4.9T, determines m*g* p = 
27r?i 2 n//iBBp=5.7. However, the relationship between 
An and B is not linear but must be deduced from the 
tilted-field experiments. We derive the actual An vs B 
curve by integrating the empirical equation of x(P) an d 
obtain An=0.69[exp(0.138£n-°- 4 ) - 1]. The curve for 




FIG. 2: Density dependence of m*g* determined by two dif- 
ferent methods. Solid data points are from tilted-field ex- 
periment with different indices: square for i=l/2, diamond 
for i=l, down triangle for i=3/2, up triangle for i= 2, circle 
for i=5/2. Parallel dashed lines indicate power law depen- 
dence of m*g* with a single exponent for all i. Their coeffi- 
cients depend linearly on i, see Fig. ^b). Open circles show 
non-monotonic density dependence of m*g* derived from full 
polarization field, B p , of in- plane field method. Inset shows 
net spin An=Pn for n=lxl0 10 cm -2 with interpolated regime 
(solid line) and extrapolated regime (dotted line). B P =4.9T 
from in-plane field method and B e:E t=6.5T from extrapolation 
of tilted-field method. A nominal m*g* ex t, slope of the dash- 
dotted line is derived from B elt for all densities and plotted as 
a thin solid line in full figure. Thick solid line represents ex- 
trapolation of m*g* to P=0 limit. Calculations from Ref. 
are shown as crossed circles. 



n=lxl0 10 cm~ 2 is also shown in the inset. The solid por- 
tion represents the interpolated regime of x while the 
dotted portion represents extrapolation of \ beyond our 
tilted field data. Requiring An=n (=lxl0 10 cm~ 2 ), we 
obtain B ext =6.5T for the full polarization field, which 
is 33% higher than the measured value of B p . This 
B ext yields a "nominal" m*g* ext = 2t:Ti 2 yl/ n B B ext =4.3, 
equivalent to the slope of the dash-dotted line in the in- 
set. 

Performing this derivation for all densities using the 
interpolation equation we arrive at an m*g* ext , which is 
plotted as a thin, curved line in Fig. It shows a quali- 
tative similarity to the m*g* p data derived from the par- 
allel field measurements, particularly the non-monotonic 
density dependence. In fact, in the low density limit 
our derived curve matches m*g* p very well, indicating 
that our empirical equation extrapolates well into this 
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regime. At higher densities considerable discrepancies 
arise, which one may, at first glance, attribute to a dis- 
continuity in An on the polarization curve, indicating the 
existence of a first order phase transition. However, given 
the good agreement between m*g* p and m*g* ext at low 
densities, and hence the absence of a phase transition in 
this regime, it is unlikely that such a transition occurs at 
higher densities. Instead, we suspect that the extrapola- 
tion of our empirical equation becomes less accurate with 
increasing density due to the increasing range of extrapo- 
lation and the solid curve of the inset bends sharply, but 
continuously towards B=4.9T. A mass increase, caused 
by the in-plane field, could provide a mechanism for such 
an accelerated bending. This effect is negligible at low 
field and hence low density, but increases rapidly with 
Hi!- Alternatively, at high polarization, terms 
of higher order in An may increasingly contribute to Xj 
leading to an accelerated bending as well. In spite of 
this discrepancy for extrapolations at high densities, the 
qualitative and partially quantitative agreement between 
m*g*ea;t and m*g* p strongly suggests that our empirical 
equation captures correctly the underlying physics of the 
system. 

Our analysis of m*g* provides a simple interpretation 
of the non-monotonic behavior of m*g* derived from the 
parallel-field data. The unusual density dependence of 
m*g* p results from a combination of the polarization and 
density dependence of \- Clearly, this m*g* does not 
reflect the actual susceptibility at any polarization value 
and consequently its density dependence cannot be used 
to assess the possibility of a ferromagnetic transition. By 
comparing measurement and extrapolation, we conclude 
that a first order transition is unlikely to occur in our 
2DES within the regime of r s studied (3-12.4). In the 
remainder of the paper, we examine x m the limit of 
vanishing polarization, which plays a central role in the 
context of a second order ferromagnetic transition. 

Our empirical formula x=m*g*=(2.73+3.9Pn)n -0 - 4 , is 
readily extrapolated to P=0 to yield the spin suscepti- 
bility, X: °f a normal Fermi liquid (see Fig. [2J). From 
n=5xl0 9 cm~ 2 to 4xl0 10 cm -2 , x = 2.73n -0 4 showing 
an enhancement factor of 1.6 to 3.6. This extrapolation 
is very reliable since it extends only slightly beyond the 
range of our data. Furthermore, from the excellent agree- 
ment at low densities between the m*g* p data and the 
extrapolated m*g* e2 , t we are confident that it provides 
very good estimates for \ t° densities as low as the MIT 
density n=2xl0 9 cm~ 2 (r s = 12.4), where x reaches about 
5.5. Whether a divergence occurs at yet lower density 
cannot be inferred from our data. 

Recent quantum Monte Carlo (QMC) calculations 
have predicted a weakly first order ferromagnetic tran- 
sition at rs=20-30, corresponding to 3.4-7.6x 10 8 cm~ 2 in 
our 2DES[6j, |7j. This density range is currently beyond 
our reach. However, our measurement of x as the 2DES 
approaches the phase transition can be compared to the- 



ory, thereby providing guidance as to the general appli- 
cability of the theoretical approach. In Fig.0 the results 
of Ref.Q are plotted as large, crossed circles i 24]. The 
general trend of increasing x with increasing r s agrees 
very well with our data although their estimates of x 
seem to be offset to slightly larger values compared to 
experiment. 

In conclusion, we have determined the polarization and 
density dependence of the enhanced spin susceptibility 
X ~ m*g* in a 2DES in GaAs/AlGaAs over a wide range 
of low densities. Our analysis provides an explanation for 
the anomalous g-factors reported recently in GaAs 2DES 
and 2DHG. The susceptibility x of the unpolarized Fermi 
liquid increases with decreasing density, in qualitative 
agreement with recent QMC calculations. 
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